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The anomalous dimensions of local single trace gauge invariant operators in M = 4 supersym- 
metric Yang-Mills theory can be computed by diagonalizing a long range integrable Hamiltonian 
& \ by means of a perturbative asymptotic Bethe ansatz. This formalism breaks down when the 
number of fields of the composite operator is smaller than the range of the Hamiltonian which 
coincides with the order in perturbation theory at study. We analyze two spin chain toy models 
which might shed some light on the physics behind these wrapping effects. One of them, the 
Hubbard model, is known to be closely related to Af = 4 SYM. In this example, we find that the 
knowledge of the effective spin chain description is insufficient to reconstruct the finite size ef- 
fects of the underlying electron theory. We compute the wrapping corrections for generic states 
and relate them to a Luscher like approach. The second toy models are long range integrable 
Hamiltonians built from the standard algebraic Bethe anstatz formalism. This construction is 
valid for any symmetry group. In particular, for non-compact groups it exhibits an interesting 
relation between wrapping interactions and transcendentality. 
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1 Introduction and discussion 



In P the one-loop spectrum of single trace local gauge invariant operators made out of the 
scalars of Af = 4 supersymmetric Yang-Mills theory was reduced to that of a nearest neighbors 
integrable Hamiltonian with 50(6) symmetry. In particular, single trace operators made out 
of two complex scalars X and Z were mapped to states in a one dimensional spin 1/2 ring, 

ti(zz ...zxz ...zxz ...zz) < — ► | tt ... TIT ... TIT ... TT) (1.1) 

Soon after it was understood that integrability persists for the full set of PSU (2, 2\4) fields 
[2j [3] and at higher orders in perturbation theory [4] where the Hamiltonian becomes long 
ranged with the range being the order in perturbation theory one considers. Later on, inspired 
by string theory data [SI El C], the full PS77(2, 2|4) Bethe equations were proposed [8] and the 
solutions to these equations are believed to yield the spectrum of generic length L operators up 
to order g 2L . At this order the interactions wrap the single trace operator and invalidate the 
use of the Bethe ansatz formalism. To achieve such remarkable point where the spectrum of 
long operators is believed to be known, a crucial step was required. Namely the idea of looking 
at operators like (11 .ip as a vacuum (the Z fields) on top of which particles (in this example the 
X fields) propagate [9]. In this language the relevant object becomes the S'-matrix scattering 
these particles, also known as magnons. This ^-matrix is S77(2|2) 2 extended symmetric [10J 
and it turns out that symmetry alone almost fixes (up to an overall function) the (4 4 ) entries of 
this matrix JTUJ [TT], 112] . The unfixed overall scalar factor has also been conjectured in [T3"l [14]. 
Knowing the S'-matrix of the theory it is then possible to write down the Bethe equations 
quantizing their momenta and, knowing the respective dispersion relation, to compute their 
energy. For example, states made out of two magnons will be given by 

£ ( e ipini+ip2n2 + S{ Pl ,p 2 )e ipin2+ip2ni ) \ni,n 2 ) + (1.2) 

where \ni, n^) represents the state with X fields in the ni'th and r^'th positions in a sea of L — 2 
Z fields. The dots correspond to a non-trivial part of the eigenstate in the boundary of the 
asymptotic region when n\ is not very far from n 2 and the magnons are strongly interacting. 
The momenta are then quantized via the Bethe equations 

e^ L = S( Pl ,p 2 ) , e* P2L = S(p 2 , Pl ) 

which physically simply mean that the phase acquired by each magnon when going around the 
ring equals the free propagation phase pL plus the phase shift due to scattering with the other 
magnons. The spectrum is then given by the sum of energies of the individual magnons as 
A — L + 2 = e OQ (pi) + €00(^2) + 0{g 2L ) where the infinite volume dispersion relation, also fixed 
by symmetry fTUJ dSl EE2] , is given by 



too{p) = y 1 + 165f 2 sin 2 - . 
The simplest possible 2 magnon state is the well known Konishi operator 

i*o Him) -1 urn (1.3) 
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Figure 1: Results and conjectures for the scaling dimension A(g) = 4+12<? 2 — 48g 4 +336.g 6 — ag 8 , of the Konishi 
operator up to four loops, when wrapping interactions first appear. The lines r\ and correspond to the two 
recent (disagreeing) computations in [flCET] for which a = 2607+28C(3) + 140C(5) and -2584+384C(3)-1440C(5) 
respectively. The first one was done using superspace techniques whereas the latter used component formalism 
making the comparison between these two laborious computations far from easy. These computations differ from 
all previous conjectures by the presence of C(5)- c± is the most recent conjecture |19j . which is based on some 
transcendentality observations and pomeron considerations and predicts a = — 5307/2 + 564£(3). Conjectures ci 
and C3 appear in |13| . The former corresponds to a — 5640 + 288£(3) and would be the value of the anomalous 
dimension of the Konishi operator if we were to believe the Bethe ansatz equation beyond their natural limit 
of validity and is therefore a very unlikely proposal [19j . C3 is the anomalous dimension whose transcendental 
part is that given by the BAE while the rational part is that predicted by the Hubbard model and has therefore 
a = 5088 + 288((3). Finally c 4 with a = 5088 would be the result predicted by the Hubbard model [5D] which 
we now know only reproduces the good BAE up to 3 loops. 



whose dimension can be computed from the known Bethe ansatz equations [HE] up to order 



At order g 8 wrapping interactions appear and the techniques at hand do not suffice to tackle 
this computation. Still there are already two possible results in the literature [lU [17] (see 
also [18]) where the g 8 coefficient was computed by direct evaluation of field theory Feynman 
diagrams. In figure [T] we plot the several computations, conjectures and speculative guesses, 
for the Konishi anomalous dimension up to four loops. 

Optimistically one might expect to find some extra integrable structure in M = 4 SYM 
which would allow one to treat the gauge invariant states beyond the perturbative asymptotic 
Bethe ansatz regime. A particularly appealing possibility would be that some extra hidden 
local degrees of freedom exist and the long range interactions we perceive would rather be the 
effect of integrating out these fundamental degrees of freedom. This scenario finds compelling 
evidence at strong coupling in [2H [221 E3J [21] , at weak coupling in [20] and for general coupling 



in [25]. 



In [2H [221 123], [21] quantum sigma models describing the S n subsector of AdS^ x S 5 type 



9 2L = 9 8 , 



A K = 4 + 12g 2 - 48g 4 + 336/ + 0{g 8 ) . 
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IIB superstring were seen to reproduce the long range conjectured AFS string Bethe equations 
[7J at strong coupling when the rapidities (0's) of the relativistic particles were integrated out 
thus leaving an effective Hamiltonian for the isospin degrees of freedom. 

In [2D] the BDS equations [4J, which are known to describe the SU(2) sector of the su- 
persymmetric gauge theory spectrum up to three loops, were shown to be equivalent to the 
Hubbard model at half filling where again integrating out the momenta (g's) of the electrons 
yields an effective long range Hamiltonian with SU(2) symmetry for the spins of the electrons. 

In the Hubbard model the effective magnons appearing in (jl.2p can be understood as bound 
states of empty sites (o) and doubly occupied sites (|). As we will describe below, if, in the spirit 
of [26l [27J [28j [29j [30] , we want to diagrammatically compute the finite size corrections to the 
effective magnon theory coming from the Hubbard model using the Luscher approach [31j [32] , 
then we need to take into account the fact that the magnon is a bound state of two fundamental 
particles rather than a fundamental excitation itself. For example, as discussed in section 12.21 
the leading finite size corrections to the magnon dispersion relation can be reproduced for any 
value of the coupling g from the expression 



where, in order to reproduce the correct result, we must use the scattering matrix between a 
magnon bound state $ and its fundamental constituents o and |. 

Curiously, Janik and Lukowski [27] computed the leading finite size correction to the Hub- 
bard magnon energy at large g = —t/U, using instead the magnon- magnon scattering matrix 
and still got a sensible result (correct up to a factor of 2 - see equation (71) in [27]). Physically 
this makes sense because at strong coupling g - which from the Hubbard point of view corre- 
sponds to weak interaction strength U compared to the electron hopping kinetic energy t - the 
magnons are a weakly bound pair of o and J and when we scatter the magnon against another 
magnon we are effectively scattering it against two fundamental particles! However, as the 
coupling decreases, the magnon-magnon Luscher formula with BDS magnon-magnon S-matrix 
will no longer provide the correct Hubbard result 0- As explained in section fl?2\ the two results 
will agree to leading order - when the correction is of order 1/g - and start disagreeing at next 
to subleading order - at order 1/g 3 . Could we be in a similar situation in M = 4 SYM? There, 
the finite size corrections [33], [311 EH] to the Giant Magnon [36] were reproduced at leading 
order J27J [28] and at next to leading order [29], [30] but no two loop computation is available. 
Bearing in mind what happens in the Hubbard model it is not completely inconceivable that at 
this loop order the Luscher results based on the lightcone S-matrix [TUJ [T2] start failing. On the 
other hand from the string worldsheet point of view this scenario would certainly be intriguing. 

In the opposite regime, at weak coupling g, expression (11.41) gives precisely the correct result 
whereas the use of the magnon-magnon BDS S-matrix is completely inappropriate because in 
this regime the elementary particles that make the magnon are highly bound. If this qualitative 

x Let us also remark that considering also the contribution from the S-matrix between magnon and bound 
state of magnons (which exist in the BDS theory) in the Luscher formalism does not cure this problem. At 
most, these contributions can reproduce the higher winding number diagrams of the fundamental constituents 
o and l only in the strong coupling regime g 1. The reason is again that in this regime a bound state of b 
magnons is almost like 2b free fundamental particles. 




) 
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structure persists in the full M = 4 theory then it explains why a naive computation of the 
Luscher terms at weak coupling seems to never yield any sort of transcendental numbers such 
as C(3) or £(5) which typically appear in the computations of [HI [T7] . The reason would be 
that in order to probe the weak coupling limit of the theory the knowledge of the magnon 
constituents would be essential 0. 

In particular it seems clear that a thermodynamic Bethe ansatz (TBA) approach to the 
BDS equations would not recover the Hubbard finite size corrections (because in particular 
the TBA approach always reduces to the Luscher formulae at large radius and as we explained 
those only work when we take into account the scattering of the magnons with the fundamental 
Hubbard electrons). In the context of the full M = 4, the TBA program is being carried out 
in [37], [38] still with inconclusive results. 

Still in the context of the Hubbard model we analyzed in sections 12.3.11 and 12.3.21 how 
wrapping interactions manifest themselves for many particle states. In the Af = 4 context this 
is an unavoidable question if we want to understand the full anomalous dimensions of small 
operators such as the Konishi operator (II. 3p . There are two different kind of effects one needs 
to take into account to study wrapping corrections to many particle states. On the one hand 
the energy of the state as a function of the magnon momenta changes when the theory is put 
in finite volume and this leads to a Luscher type correction which for the Hubbard magnons 
reads 

1 M . , M 

6E Luscher = -j2Yl / ^(^9)-4W)e i(9 - Wi n4:>-?) +c - c -- 

n=l a =o,\ JCn 71 m=l 

generalizing (11.41) . Led by the striking simplicity of this expression we conjecture a generaliza- 
tion of the Luscher formula for many particle states in integrable two dimensional models in 



The second effect we need to take into account is the fact that due to the wrapping in- 
teractions the quantization conditions - that is the Bethe ansatz equations - for the magnon 
momenta receive corrections and thus the momenta are slightly shifted when wrapping inter- 
actions are taken into account. For example, BDS equations [I] can be dressed to 

+ \ L — j_ • 



n 



e 



X„ J U n Um % 

in such a way that the leading wrapping interactions are taken into account. Here cj) nm is a 
wrapping dressing kernel described in section 12.3.11 

As explained in section 12.3.31 a particularly curious feature of the computations of the 
weak coupling finite size corrections is that the leading wrapping correction, of order g 2L , to 
a state whose magnons' momenta are pj ~ p^p + g 2 Pj + • ■ ■ + g 2L pf J ^ only depends on the 



2 On the other hand, if the magnons in the light-cone gauged string theory can be thought of as fundamental 
particles - contrary to what happens in the Hubbard model - then a priori we should indeed sum over all these 
infinite number of bound states and of course infinite sums can eventually lead to transcendental numbers. This 
seems the only way out to find the good transcendental Konishi anomalous dimension from a weak coupling 
computation based on the Luscher formulae. We acknowledge N. Gromov, V. Kazakov and K. Zarembo for 
discussions on these issues. 
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leading values . This is natural from the point of view of the Luscher computations which 
are basically a smart way to organize the two dimensional perturbative expansion using the 
two dimensional S'-matrix. In this formalism the exponentials which appear in the several 
integrands are automatically of order g 2L and thus the rest of the integrand, including the 
S-matrix and the dispersion relations, can be treated at g° order. In particular we can easily 
compute the wrapping correction to any many magnon BDS state without performing any 
iteration of the BDS equations up to order g 2L . In M = 4, this should be related to the fact 
that to compute the wrapping corrections to the Konishi operator we isolate the appropriate 
wrapping diagrams and only keep some information about the lower order graphs (T6J Ej- If 
the four dimensional wrapping Feynman diagrams in [HJ [TTJ, [18] could be re-written in a two 
dimensional language this could be the key to understanding the structure of an hypothetical 
hidden level of fundamental particles. 

In section [3] we consider a completely different type of (toy) model where wrapping inter- 
actions are under control. Namely we study long ranged Hamiltonians coming from a transfer 
matrix algebraic construction a la Leningrad school. In the algebraic Bethe ansatz formalism 
the fundamental object is the transfer matrix, which is a trace of product of i?-matrices, 

T(u) =tT aux R L (u)®---®R 1 (u), (1.5) 

where the /^-matrices are (simple) matrices obeying the Yang-Baxter relation and acting on 
the product of an auxiliary space (common for all i?-matrices in this expression) and a physical 
space h n . The full Hilbert space of a L-site spin chain is given by the tensor product 

H = h l ®---®h L . (1.6) 

The transfer matrix is then an operator acting on the full Hilbert space (and by definition 
scalar w.r.t. the auxiliary space). The algebraic Bethe ansatz program yields us the spectrum 
of the transfer matrix T(u). The diagonalization of this object is of great importance because, 
as we review in section [31 by taking n derivatives of the logarithm of the transfer matrix T at 
n = Owe can generate integrable Hamiltonians of range n. In particular if we take more than L 
derivatives of this object we will start generating Hamiltonians which are long ranged, contain 
wrapping interactions and still, by construction, are integrable and diagonalized through a set 
of Bethe equations. For example if we consider 



1 00 „2n Jn 

n=l 



+ h.c. (1.7) 

A=0 



then we will have an Hamiltonian which at order g 2n is of range n and whose spectrum is given 
by a sum of individual dispersion relations plus a wrapping correction which starts precisely at 
order g 2L , 

M 

wrapping (pi, . . . ,p M ) 

0=1 

This behavior is probably completely generic and we considered Hamiltonians of the form (jl.7p 
with SU(2), SU(N) and SL{2) symmetry. 
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In the SL(2) case we found the following curious behavior: suppose we consider a Hamil- 
tonian of type (jl.7p with a n some algebraic numbers. Then the dispersion relation truncated 
at a given order g 2n is a rational function of these algebraic numbers and of the Bethe roots 
Uj = ~ cot y , which are quantized via a set of polynomial equations. Thus the contribution to 
the spectrum of the X]j=i e (p?) t erm w iU be given by some algebraic quantity. On the other 
hand, precisely at order g 2L the wrapping corrections enter the game and those are given by 
an infinite sum (13.171) of algebraic functions of the Bethe roots Uj. Typically they will give rise 
to transcendental contributions! 

As an example, in the table below we listed a couple of energies of some two magnon states 
up to order g 2L and for the simplest choice of Hamiltonian with a n = 1 we get 



L 


«1,2 


e(pi) + e{p 2 ) 


E wra pping (pi j P2 ) 


3 
4 

6 


5 10 
1 _L V7 
3 6 

±1 + 

3= 2 T 2 


llg 2 , 13V3S 4 5g 6 
8 ~ 32 6 
5g 2 g 4 19g e > , 3<? 8 
4 8 24 ' 2 
a 2 + 9 4 9 6 A /o„8 4 9 10 5V2g 12 
y"^V2 3 VZ i/ 5 3 


i(_26 + 37r 2 -4C(3))/ 

/ 41 5tt 2 tt 4 \ 8 
^ 32 48 360 J y 

^(7-4C(3)-2C(5))^ 12 



(1.8) 



It would be very interesting to explore this connection between transfer matrices of non-compact 
groups and transcendentality in the context of M = 4 supersymmetric Yang-Mills theory. 
Perhaps this could provide us with important hints about the origin of the dressing factor 
which is populated by transcendental numbers. 

This paper is organized as follows: After introducing the Hubbard model in section [2] and 
reviewing the magnon description in section [2~T1 we re-derive the exact finite volume dispersion 
relation by means of Feynman diagrams in section I2.1.11 In section 12.21 we explain how the 
leading finite size corrections can be understood from a Luscher type approach. In section 
12.31 we study wrapping effects for general many particle states. In particular we review how 
the BDS equations follow from the Lieb-Wu equations, explain how then can be upgraded to 
include the leading wrapping corrections (section [2.3. ip and analyze the analogue of the Luscher 
formulas for many magnon states (section [2. 3. 2\i . Section [3] is devoted to the study of integrable 
long ranged Hamiltonians derived from an algebraic Bethe ansatz formalism and in section 13.11 
we explore some generalizations of this construction and present a curious non-compact long 
ranged Hamiltonian where transcendentality and wrapping are intimately related. 

2 The Hubbard Model 

The one dimensional Hubbard model describes spin 1/2 electrons moving in a periodic lattice 
with L sites. The electrons can hop between neighboring sites and there is a repulsive (or 
attractive depending on the sign) potential when two electrons (with opposite spin) occupy 
the same lattice site. Obviously, due to Pauli exclusion principle no two equal spin electrons 
can ever occupy the same position. At half filling, when the number of electrons equals the 
number of sites, each electron will tend to occupy a site in the lattice due to the repulsive 
potential. We can then study an effective Hamiltonian for the spins alone [39]. It will be a long 
ranged Hamiltonian where the interactions correspond to virtual processes where electrons hop 
there and back eventually changing spin in the process. In [20] this effective Hamiltonian was 



S 



identified with the long range Hamiltonian of Af = 4 supersymmetric Yang-Mills theory. This 
identification is correct up to three loops but fails beyond that. Still, this is an instructive toy 
model since wrapping interactions, due to electrons making loops around the ring, are perfectly 
under. In this section we will study them, give them a diagrammatic description and understand 
how they fit into the usual field theoretical Luscher treatment of finite size corrections. We will 
also understand how to modify the effective Bethe equations for the magnons of the effective 
spin theory so that they reproduce the (leading) wrapping effects. 

A quite useful alternative description of the relevant degrees of freedom in the Hubbard 
model is obtained by performing a so called Shiba duality. It amounts to thinking of the Hilbert 
space as that where N — M vacancies o and M double occupancies | move in a ferromagnetic 
vacuum with L up spins^). In this description the Hamiltonian reads 

L L 
i=l a=o,l i=l 

where, following [20], we have introduced some extra twists <p a in the Hamiltonian which can be 
thought of as a sort of magnetic flux inducing additional phases in the electron wave function 
as it moves around the chain. As explained in [20] and reviewed below these twists can be used 
to delay the wrapping corrections to the effective spin theory. 

The complete spectrum of this Hamiltonian can be obtained as 

N 

E = J2<Qn) , e(q) = -2tcos(q) (2.2) 

n=l 

where the momenta are quantized through the solution of the Lieb-Wu [10] equations, 

J(q n —(f>o)L 



e 



n' [ Ui — 2osin(o r) ) — ill 
Uj - 2g sm(q n ) +i/2 

e »(^i)LTT «j-2gsin(g B )-V2 = -pj Uj — Uk — i j = l,...,M,. (2.4) 

„ = i u i ~ 2 9 sm(g„) + z/2 ^ u 5 - u k + % 

where 

g = -u- 

In section 12.31 we will review [20J how, eliminating the electron momenta q n from these 
(twisted) Lieb-Wu equations, we obtain an effective set of (twisted) BAE for the spin degrees 
of freedom u which are precisely the (twisted) BDS equations [1]. In that section we will analyze 
wrapping corrections in full generality. We will see for example, that perturbatively in g the 
effective (twisted) Bethe ansatz equations are valid up to order g 2L for the choice of twists 



y = mod ^ (2.5) 



3 The number of vacancies o and double occupancies | is related to the number of up and down spins as 
7V| = N u L - N = JV T . 
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and the g correction to any perturbative state will be given solely in terms of the position 
of the Bethe roots Uj to leading (g°) order only. This is probably a simplifying feature of the 
weak coupling finite size corrections which is likely to be present in the Af = 4 supersymmetric 
spin chain. 

We will also study the generic case where (12.51) does not hold since it will provide us with 
a nice toy model to understand Luscher type corrections for many particle states. 

In the following two sections we will consider a much simpler setup which however captures 
most of the relevant physical information. Namely we will consider the very simple configuration 
with a single vacancy o and a single double occupancy J, that is N = 2 and M = 1. 

2.1 The magnon 

A magnon in the Heisenberg XXX spin-| chain, 

H xxx = - ^2 (1 — ' <?n+l) (2-6) 
n 

is the lowest lying excitation above the ferromagnetic ground state. It is a plane wave state 

E eiPV «it---t) 

n 

of one down spin in a chain of up spins, with excitation energy 

e(p) = 1 — cosp . (2.7) 

In the Hubbard model such plane wave is not an eigenstate of the Hamiltonian but there is a 
close analogue of this state when the empty site excitation o and the double occupancy J form 
a bound state (note that a o and a J in the same site is precisely the same as a spin down). 
More precisely, as reviewed in appendix A, we can consider the following half-filling state 

|*) = 5>(n,n') (al x al >0 ) | t ■ ■ ■ t) ■ 

n,n' 

with ip(n, n') being a superposition of plane waves with momenta q and q' . By acting with the 
Hamiltonian on this state we can see that there exist bound state solutions witb@ q,q' = | ± i(3 
where p is the bound state momentum while (3 dictates the exponential damping of the wave 
function away from n = n'. The energy of such states, which we also call magnons, equals 

e(p) = -4t cos | cosh /? . (2.8) 

In infinite volume we find (3 = f3oo{p) where 

4( ? cos|sinh/5 00 (p) = -1 (2.9) 

4 Throughout this paper we always use this definition of p which seems the most natural one from the 
Hubbard point of view. To make contact with the standard notations in the AdS/CFT literature we should use 

Phere — Pusual "T 
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so that e(p) = eooip) with 



Coo (p) = -JlP + 16t 2 cos 2 | = -£/ (1 + V(cosp + 1) + . . .) (2.10) 

which at weak coupling g has the same — cos p dependence as fl2.7jl . This is expected from 
the known result that perturbatively in small g = —t/U the Hubbard model at half filling is a 
long-ranged Hamiltonian whose leading term is the Heisenberg spin chain. 

The magnon state can also be described by the triplet (q,q*,u) satisfying the Lieb-Wu 
equations 



e i{ q -^o)L = u-2gsm(q) -i/2 e i{ q *-<t> )L = u - 2gsm(q*) - i/2 

u — 2g sin(g) + i/2 ' u — 2g sin(g*) + i/2 ' 

c i(^~ML u ~ 2 9 sin(g) i/2 u 2g sm(q*) - i/2 = ^ 
u — 2g sin(g) + i/2 u — 2g sin(g*) + i/2 

In infinite volume, the l.h.s of the first two equations is and oo for a bound state with complex 
momentum q — | + i(3. This fixes 

u = 2gsinq + ^. (2.11) 

The reality of u implies 

p 

u = 2g sin- cosh Poo(p) = u^p) (2.12) 

and the condition ( 12. 9ft which gives the dispersion relation ( 12.101) . 

When the state is put in finite volume, equation ( 12.81) is still valid but the expression ( 12. 9ft 
for P(p) gets modified to 

a p ■ v. a sinh/3L 

Ag cos - smh p = j— (2.13) 

2 cosh f3L — cos (jp — 2<j) ) 

In appendix A we derive this equation from the direct study of the magnon wave function in 
finite volume. Naturally, the same result can also be obtained from the Lieb-Wu equations 
( 12.41) . Indeed these Bethe equations are exact for any L since the interactions of the elemen- 
tary particles are ultra local. The leading finite size correction to the magnon energy is then 
exponentially suppressed in the ratio of the system size L by the the physical size 1/fioo of the 
bound state, 

e(p) - eoc(p) = Ue- Mp)L 2t&nhf3 O0 (p) cos - (p - 2<\> ) (2.14) 

For the particular choice of twists (12.51) . such that ^°~^V L = — 1 ; this leading correction 
vanishes, and instead we have 

e(p) - eoc(p) = t/e- 2/3o °^ L 2tanh/? 00 (p) (2.15) 

In the next section we will recover these known results from a diagrammatic point of view. 
This will turn out to be very useful to understand how to recover the finite size corrections 
from the effective theory point of view, that is, from the BDS language. In particular we 
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( w >?) 

\ = l 

7 u— e(q)+ie 
(w,g) 

„N. = J 

' o>— £(g)-Ne 

Figure 2: Feynman rules for diagrammatic computations in the Hubbard model. Each ele- 
mentary particle (o and |) has a non-relativistic free propagator (solid and dashed lines) and 
there is only a quartic interaction vertex. Loops carry an extra minus sign as the elementary 
particles are fermions. 




will understand that if we were given solely the BDS Bethe equations we would not be able 
to recover the finite size corrections using a Luscher type approach [2H E2] except at strong 
coupling. Instead, the knowledge of the magnon constituents when computing the virtual 
processes wrapping the space-time cylinder will turn out to be essential. The diagrammatic 
language seems therefore promising to try to learn some lessons about what to expect for 
the M = 4 SYM chain if this chain most fundamental description comprises extra degrees of 
freedom [2H M, E21 M, M, EE] . 

2.1.1 Diagrammatically 

In this section we shall explore the field theoretic description of the Hubbard model defined by 
the Hamiltonian (12.1 1) . This will allow us to use the powerful diagrammatic techniques of field 
theory to obtain the finite size effects from loop diagrams with topological winding around the 
compact direction. We start by writing the action of the theory 

S= I dtY,{£o + C\+Cint) 
J n 

where the free part is given by 

At = - {a* np d t a n ^ - a nya d t a* n cT ) + t (e^ a* n ^a n+ x, a + e" 40<T a* j(T a n _i i(T ) 
and the interaction by 

The elementary excitations o and J have a non-relativistic propagator 



{Ta n ^(t)a* n , t(T (t') ) 

where 



2tt J_ n 2n 10 — e(q) + ie 



e(q) = — 2t cos(g) . 



The interaction term in the action gives rise to a quartic coupling with coupling constant U. 
The Feynman rules are summarized in Figure [2j In a non-relativistic theory the number of 
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Figure 3: Feynman diagrams for the two point function of the composite operator \. 



particles is conserved and this greatly simplifies the field theoretic perturbative expansion of 
the theory H2J H3]. Diagrammatically, this fact stems from the retarded nature of the 
propagators which gives them an orientation. 

In order to find the two particle spectrum we consider the two point function 

(TXn(t)Xn>(t>)) = J |1 J"j£e^*™^G(n,p) 

of the composite operator 

Xn(t) = a n>0 (t)a nt i(t) + a* no (t)a* ntl (t) 

In particular, to find out possible bound states (magnons) we should look for poles of G(Q,p) 
thus obtaining the dispersion relation Q(p). Notice that there is a big arbitrariness in the choice 
of the composite operator x- The only requirement is that the state it creates has some overlap 
with the magnon wave function. 

The propagator G(Q,p) can be computed diagrammatically. It is given by the sum of 
all diagrams describing the two elementary particles moving freely and interacting k times as 
shown in Figure [3j More precisely, it is given by 

oo 

G (Q,p) = [Go (V,p)] k+1 Ht/) fe = . r _ 1( * : (2.16) 
where the free propagator 

can be computed by residues, 

G (n,p) = 1 = 

K ,P) ^ - 16t 2 cos 2 | 

Thus, the full propagator reads 

G(0,p) 



2tt J _ n 2tt lu — e(q) + it Vt — oo — e(p — q) + it 



^9? - 16t 2 cos 2 \-U 
and the magnon bound state corresponds to the pole at 



O = + 16^ 2 cos 2 | 



13 



which is precisely the infinite volume result (I2.10p . 

With periodic boundary conditions of size L the free propagator is changed. The natural 
way to account for this effect is to sum over all possible windings of the loops in Figure [3] 
around the compact circle. The winding numbers are a topological property of the Feynman 
graph that can be assigned to any of the propagators forming each loop. We shall compute the 
graph assigning the winding number m always to the particle o. Furthermore, equation fl 2 . 1 6 [) 
remains valid provided 



G (n,p) = - Yl J 2tt 7 - e(q) + ie ft - a - e(p - q) + ie 

m=—oo v ' v ' 



e im(q-4> )L 



^/Q 2 - lQt 2 cos 
where, if we parametrize Q as 



2 P 
2 



ft = -At cos ^ cosh/?, (2.17) 
we get 

The full propagator then becomes 

g(n, P ) - iF ^ p) 

v ' y/n 2 - m' 2 cos 2 1 - uf l (n,p) 

and the magnon pole sits at ( 12.17P with /3 determined from the equation 



At cos - sinh (3 = UF L (ft, p) 



which is precisely (12.131) . 



2.2 Luscher in Hubbard 

Luscher developed a general formalism to determine the leading finite size corrections 

in quantum field theory. In particular, he studied the change in the dispersion relation of 
one particle when one imposes periodic boundary conditions. Remarkably, he found that the 
leading correction to the energy e(p) of a particle with momentum p was completely fixed by 
the particles infinite volume dispersions relations and S-matrix. The idea is that the on-shell 
dispersion relation is defined by the pole of the propagator. Following the standard notation 
of relativistic field theory we can write 

iG- 1 (ft,p) = -ft 2 + 4(p) + S (ft,p) = (2.18) 

where 

S(ft,p) = E 00 (ft,p) + 5E (ft,p) 
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is the particle's self-energy, whose infinite volume part vanishes with zero derivative on the 
mass-shell Q = e^p). From (12.181) with Q = ^{p) + Se(p) we obtain 

^-£ + (£)(© + (^-0i(S ,+ - (219) 

with all the quantities computed at the infinite volume mass-shell Q = e 00 (p). The self-energy 
correction can then be related to the S-matrix. This is achieved by evaluating the self- 
energy diagrams with winding around the compact direction by deforming the integration over 
the loop momenta to pick the on-shell pole of the wound internal propagator and obtain an 
integral of the forward scattering S-matrix [3TJ E2J EZj- More precisely, the leading finite size 
correction to e a (p) is given b}@ 

\ I ^^D-^K^-^'W) (l-S a a b b (p,q))+c.c. 
Jc b 

where S^l(p,q) is the S-matrix for forward scattering of particle a with particle b, F b = ±1 
encodes the bosonic/fermionic nature of particle b and we consider only the contribution from 
diagrams with winding number ±1. For usual relativistic theories the contour C is given by an 
integral over the possible S'-matrix poles plus an integral over the real axis. The latter describes 
a quantum loop and is absent in our case where the underlying theory is non-relativistic [HI 
H2J H3] . The former can be simply computed by residues and is denoted by /z-term. Let us now 
apply this general formalism to the Hubbard model and find the leading finite size correction 
to the magnon dispersion relation. As we saw in the previous sections the magnon is a bound 
state with a finite size 1//3qo- Therefore, its energy in a finite system has a leading correction 
of order e _/3ooi , except for a particular choice of twists where this correction can be delayed to 
order e~ 2/3ooL . The diagram leading to the first wrapping correction corresponds to the splitting 
of the magnon into its fundamental constituents each one going around the space-time cylinder 
in opposite directions and meeting on the other side of the cylinder as depicted in figure HI The 
Luscher /i-term then reads 

Hp) = E (I Res « elLiq ~ M ^ - 4(p)) 4'>' 9) + c - c ) ( 2 - 2 °) 

where q and e(q) are the momentum and energy of the elementary particle o going around the 
loop in figure HI (p) is the infinite volume dispersion relation of the magnon and the residue 
is taken at the pole of the S'-matrix. The S-matrix between magnon and elementary particle 

5 In [311 132j this formula was derived for relativistic theories and p = 0. This amounts to computing the 
corrections to the particle's mass. In the process of derivation the propagator is wound around the spacial 
circle originating a factor of cos(qL) multiplying the infinite volume propagator. Since the spacial momentum 
p vanished the problem was isotropic with respect to the spacial directions and this factor could be simply 
replaced by 2e lqL . In the case p ^ we must treat each exponential separately which amounts to summing the 
contributions from the virtual particles going parallel and anti-parallel to the physical particle. In most cases, 
including those considered in this paper, these symmetrizations will simply lead to computing the real part of 
the result obtained keeping one exponential. 
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Figure 4: The leading finite size correction to the magnon energy is given by the Feynman 
diagram where its two elementary constituents split and merge after winding the spacetime 
cylinder once. After cutting the wound loop, putting one elementary particle on shell, this 
diagram gives rise to the Luscher /z-term. 



(o or |) can be read of from the Bethe equations (12.41) 



clo ( \ _ s U (n a) _ u 00 {p)-2g S m(q)-i/2 



with 

The pole condition 
has the simple solution 



WooO) = 25(sin - cosh/3oo(p) . 



25(sin(g) = ^(p) + i/2 



V 



Putting everything together we recover the result (12.141) . For the choice of twists (12.51) this 
leading term vanishes and one needs to consider the contribution coming from diagrams with 
winding number ±2 [§, 

6e(p) = Res, j 2L ^ (e'(q) - 4(p)) Sft(p, q) + c.c. 

which gives (12.151) . 



6 For generic twists, the subleading correction (of order e _2 ^ ooL ) to the magnon dispersion relation depends 
on all terms in (|2.19p . including the second derivative of the infinite volume self-energy which is not a physical 
on-shell observable. However if SH vanishes at order e _,9ooi then we only need to compute the first term, that 
is, the correction to the self energy due to winding 2. 
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2.2.1 Luscher with magnon-magnon S-matrix 

In the previous sections we saw, by many means including a Luscher type computation, that 
the leading finite size correction to the magnon dispersion relation is given by [20] 

g e (p) = Ue~ 2l3 ™ L 2tanh[3 O0 

when the Hubbard twists are chosen as in (12.51) . This expression is valid for all values of the 
coupling g. In particular, at strong coupling we find [34] 

f U p U , p \ L p ( 1 2 p \ 

oe(p) = sec — I —- sec — h . . . exp sec - 1 sec — h . . . . (2.21) 

\ 2g 2 64c/ 3 2 J F 2g 2\ 96g 2 J v ' 

In [27] - in the process of computing the leading //-term prefactor for the AdS§ x S 5 giant 
magnon - Janik and Lukowski also determined what the contribution would be for the BDS 
scenario. They found 

x f \ U p i 
oelp) = sec -e 2 » 2 + . . . 

9 2 

where we adapted their results to our normalizations. This result was obtained applying Luscher 
formula with the magnon-magnon S-matrix and thus describes a different physical process from 
the point of view of the elementary particles. Surprisingly, the two results agree up to an overall 
factor of 2. Physically one can try to justify this result using the fact that for large g the magnon 
is a weakly bound state of the elementary particles o and J. Thus, it is not unreasonable that 
the effect of a magnon loop around the spacetime cylinder is almost the same as a loop of 
the elementary particles with winding number 2. In any case, the accidental nature of this 
agreement is confirmed by its limitation to the strong coupling regime. In appendix C we 
repeat the computation of [27J, expanding the result further in 1/g, 

x t \ sof U P , iU + P 2 P , \ L P (, i P P , \ 
oeip) = yi sec — tan - sec — h . . • exp sec - 1 tan - sec — h • • • . 

K J V 9 2 g* 2 2 J P 2g 2 \ 2g 2 2 J 

(2.22) 

So that the result starts differing (apart from the factor of 2) at the next order in 1/g. If we 
fix L/g and expand both the prefactor and exponent then we can easily see that the 1/g 2 term 
drops out and the mismatch is delayed to the next to leading order in 1/g, as mentioned in the 
introduction. This limit is the one usually considered in the AdS/CFT context. 

The BDS magnon effective theory, reviewed in the next section, contains bound states of 
magnons [B] which are simple generalizations of the Bethe strings in the Heisenberg model. 
In [28] the Luscher term accounting for the finite size correction to the (dyonic) Giant magnon 
[44J dispersion relation was written including the virtual exchange of all these bound states of 
magnons. Thus, one might question if the sum of the contributions of all these bound states 
will correct the result (I2.22p to give the exact result (12.211) . This is not the case. In fact all 
diagrams present in the theory are those in figure [3] and they do not, in any sense, describe 
virtual exchanges of magnons and magnon bound states. At most, at very large g the loops 
of b magnons could be mimicking 2b loops of fundamental particles and thus, the best we can 
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expect from summing all the possible magnons is to recover the full finite volume result but 
always for g — > oo. 

Furthermore, in the weak coupling regime the magnon-magnon computation gives the right 
g 2L coupling dependence but misses completely the momentum dependence of Se. 

This result puts in question the validity of direct application of Luscher formulas to the 
giant magnon of M = 4 SYM. In particular, if the giant magnon is not an elementary particle 
(like the 2 electrons bound state in the Hubbard model) of the worldsheet theory then its finite 
volume energy will be sensitive to the elementary excitations of the theory and can not be 
recovered just from the magnon-magnon S-matrix in infinite volume. On the positive side, this 
makes the wrapping effects a window into the most elementary level of the theory. 



2.3 Magnon effective Bethe Ansatz equations and generic wrappings 



In this section we will review the results of [20] and see how the Bethe equations for the effective 
system of spins at half-filling coincides with the (twisted) BDS equations. We will then study 
the wrapping interactions for a general state. In particular we will understand how to dress 
the BDS equations in such a way that they include the wrapping interactions. We will also see 
that the leading order correction to the energy of a generic M magnon state can be given a 
Luscher type diagrammatic interpretation. 

At half filling the Bethe ansatz equations ( 12. Ah are given by 



e 



i(<ln-<t>o)L 



' l uj- 2gsin(q n ) -i/2 



l u j -2gsm(q n )+^/2 , 1 ' ' 



M 



e 



i{ln+M-<t>o)L 



TT uj - 2g sm(q n+M ) - i/2 
}=i u i ~ 2 9 sin(g n+ Af) + i/2 ' 

C i(^-ML TT u i ~ 2gsin(g n ) - i/2 Uj - 2gsm(q M+n ) - i/2 _ -pr Uj - u k - i 
_i u j - 2 9 sin(g„) + i/2 Uj - 2g sm(q M +n) +i/2 £-1 u . _ Uk _|_ j ' 

where we explicitly split the 2M momenta q n into two equal groups so that all indices range 
from 1 to M. A state with M magnons is a state where the 3M Bethe roots organize in triplets 
of Takahashi states 

(q n , q n+M , u n ) , (2.26) 

with 

_ * 

Qn+M — Q n 5 

while the real u root is given by 

i i 
u n ~2gsmq n + - = 2gsinq* n --. (2.27) 

The last condition is necessary once we allow q n to have a positive imaginary part (and thus 
q n +M = Q n to have a negative imaginary part). In this case, the l.h.s of equations (12. 23ft and 
(12.241) . respectively, vanishes and diverges exponentially with the system size L and condition 
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(12.271) is required to obtain the same behavior for the r.h.s. It is also clear - and shown in 
[20J - that if we multiply the equations for each element of the Takahashi triplet we obtain the 
twisted BDS equations 



m _ . 

( . ... U n u m ~T 1 



J( Pn -<po-cf> t )L ^ TT "n ^ (2.28) 

parametrizing the momenta as 

Pn . ■ n 
q n = — + iPn ■ 

The dispersion relation e(q n ) + e(g n+ A/) as function of p and the relation between p and the 
Bethe roots u is to leading order exactly as before, as it follows simply from computing the real 
and imaginary part of (12.271) . namely [3 n ~ /3oo(Pn) given in ( 12. 91) . 



u n ~ Moo(Pn) = 77 tan^Wl + 16g 2 cos 2 ^ , (2.29) 



and E ~ Y, e oo(p n ) with 



eoo (p) = - A /f/ 2 + 16t 2 cos 2 ^. (2.30) 



2.3.1 Corrected BAE 

In this section we will study expressions (12. 28|2. 29112. 301) in greater detail. That is we will 
understand how these relations get modified once the leading finite size effects are taken into 
account. To do so we will study the leading wrapping effects for a generic many magnon state. 

Physically there are two sources of corrections to the energy of a many magnon state. On 
the one hand, the energy of the state as a function of the magnon momenta changes when the 
state is put in finite volume. This will lead to a Luscher type correction. On the other hand, 
the periodicity condition, that is the BAE, are corrected and thus the quantized momenta 
are shifted. Diagrammatically this last effect would be due to new wrapping virtual processes 
correcting the magnon S-matrix rather than to the usual self energy virtual graphs present for 
the Luscher type contribution. 

As for the single magnon case we will see that the bound state structure of the magnon 
must be taken into account to reproduce the proper finite volume results. 

Needless to say, instead of correcting the effective BDS equations we could simply use the 
exact Lieb-Wu nested Bethe equations! The point is that we want to understand how the 
corrections to the Bethe equations of effective spin theories look like. This might be useful 
if, as discussed in the introduction, the M = 4 SYM long-ranged Hamiltonian stems from an 
underlying Hubbard like description. 

We want to eliminate the magnon momenta q n from the Lieb-Wu equations thus obtaining 
an effective equation for the magnon rapidities u n . As we saw, to leading order we simply have 
(12.271) but since we want also to keep track of the leading wrapping corrections to the effective 
equations we should instead write q n = ^ + i(3 n and 

u n = 2g sin q n + - + A n = 2g sin q* n - - + A* (2.31) 
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where A n is a small quantity which can be computed from the equations (12.231) and (12.241) for 
the magnon momenta. This is done in appendix B. Since u n is real we can compute both u n 
and (3 n from the knowledge of the real and imaginary of the small quantity A n . In particular, 
taking the imaginary part of (12.311) . we obtain that Sj3 n = (3 n — /3oo(p n ) is given by 

5(3 n ~ 2tanh/3 00 (p n ) 3 (A n ) (2.32) 

where as before Poo(p) is defined through (12. 9j) . To proceed we introduce the notation = 
+ if3oo(Pn) and = 2g sin — | so that 



u, 



u™ + 2gi cos q™ 5f3 n + A n + 0{e~ 2(3L ) , (2.33) 



Next, as explained in greater detail in appendix B, we multiply the three Lieb-Wu Bethe 
equations for the Takahashi triplet (I2.26P and expand using ( 12. 31 [) and (12.331) to find 

u~ - u™ + i 



FT ^ -^—^ e l<t> ™ + 0(e~ 2/3L ) , (2.34) 

J. J. 7/OO 7/ 00 , \ / 1 \ I 



m^=n n "' 



where the wrapping dressing kernel 4> nm is given b}0 

3(A n ) (I 2 p n 



. — tan 2 — + - (n «-> m) (2.35) 

u n a m ) "T 1 \ a n * a n a m / 

The dressed Bethe equations (12.341) can trivially be solved perturbatively provided a solution 
p°^ to the original BDS equations is given. In this case which was a function of p n , can be 
expanded around the value p n = p^. The value of the Bethe roots for these values of momenta 
are denoted by u^ 3 . We stress again, the u^ 3 are the values of the Bethe roots obtained via 
the BDS equations whereas are functions of a free variable p n . Writing = + 5u n we 
easily see that the leading order shifts to the Bethe roots due to the inclusion of the dressing 
Kernel reduces to the simple linear problem 



L8p n - 



2 (8u m - 5u n 



^ (u~ - u^) 2 + 1 t-f 



with 5u n = ( ^r 2- ) $Pn- Having found 5p n and 53 n we want to compute the shift to the 

\ ap n I _ oo 
^ ' Pn—p n 

many particle state energy 

SE = J2 ~ At cos V cosh ^ n + At sin \ cosh ) 

n 

due to the wrapping effects. As mentioned above, this expression is non-zero due to two 
completely distinct type of corrections. 

On the one hand, as we just saw, the momenta are quantized differently and thus we will 
have a contribution due to the displacement of the Bethe roots when wrapping interactions 



7 If the twists are chosen as in (|2.5[) the imaginary part of A„ becomes of order e 2f3L and we need to expand 
further. In this case we obtain (|B.6|) in Appendix B and (|2.33l2.34p hold to order e~ 3/3L . 
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are taken into account. For many particles states like the ones we are now considering these 
corrections must be taken into account. 

On the other hand the functional dependence of the energy on the momenta {p n } changes 
when we put the system in finite volume. More precisely (3 n will be given by the infinite value 
expression f3oc(pn) plus the correction 5f3 n which will in general depend on all the magnon 
momenta in an entangled way. This contribution is precisely the analogue of the Luscher 
corrections described in the previous sections for the single magnon case. We will discuss these 
corrections in greater detail in the next subsection. 

Putting these two corrections together we immediately get 

SE = (~U5Pn + ^j^-Spn) ■ (2.36) 

Roughly speaking we could say that the first term is of Luscher type and accounts for virtual 
processes correcting the magnon dispersion relations while the second term stems from correct- 
ing the magnon S-matrix and therefore the BAE. It would be interesting to provide <p nm with a 
diagrammatic interpretation. Moreover, in the thermodynamical Bethe ansatz approach to the 
computation of finite size effects in relativistic integrable models, renormalized Bethe equations 
for the positions of extra zeros in the TBA F-system appear If the same structure emerges 
for the AdS/CFT TBA equations then the form of the above dressed equations (I2.34p might 
provide some hints about the possible aspect of such dressed equations. 

2.3.2 Meeting (and generalizing) Luscher 

In this section we want to analyze the first correction SE Luscher = —U Y^ n =i ^° ^ ne ener gy °f 
a M magnon state when put in finite volume and provide it with a simple physical diagrammatic 
meaning^. Using (12.321) and the expression for A n in appendix B we find 

TT M M oo _ oo I ' 

6E Luscher 4V2 tanh/UPn) e<*«*-M-+-) L TT ^ * % + ^ 
which can be written as 

M M 



5E 



Luscher 



\ E E / ^ - e 'oo(^)) el(q ~ ML II S l>rn, q) + cc. (2.37) 

n=l a=o,\ n m=l 



where C n encircles the pole of Sf°(p n , q) at q — + if3oo(Pn) in the counter-clockwise direction. 
Obviously this expression resembles the single magnon Luscher formula (I2.20p used before and 
can be thought of as its many particle generalization. Physically it represents the correction to 
the state self-energy due to the process where a virtual particle with momenta q goes around 
the cylinder scattering with all other physical excitations as depicted in figure [51 

Due to its strikingly simple form one might try to guess what the many particle Luscher 
formula for generic quantum integrable two dimensional field theories with factorized scattering 



3 This section, together with appendix D benefited largely from discussions with K. Zarembo. 
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Pi P2 



Figure 5: The many particle state is corrected due to interactions with a virtual particle going 
around the spacetime cylinder. For integrable theories the correction to the energy of the state 
is expressed in terms of a product of factorized scattering matrices between the virtual particle 
and the various physical particles. For non-diagonal scattering this product defines a transfer 
matrix, a central object in quantum integrability. 



could be. A likely candidate for such expression for a state with M particles with polarizations 
ax, . . . , Qm and momenta pi, . . . , pu would be 



M 



5E Luscher = E f^i( £ 'M- £, a n (Pn)) e<*-+*) L (-l)^ (2.38) 

n=l {b u ...,b M } J 

s££(pi, q)... s^iiM q ) - < 2 . . . <t 

where we sum over the fundamental polarization bj but in principle allow the physical particle to 
be bound states in which case some of the indices a n will be bound states indices as in (I2.37p . In 
this case, the corresponding S-matrices in this expression should be the usual fused ^-matrices. 
For fermionic virtual particle we include the standard —1 factor from the loop and in case 
magnetic fields are coupled to any of the particles a corresponding twist is included. Notice that 
we allow the polarization of the virtual particle to change as it scatters with each of the physical 
particles. This is not in contradiction with (12.371) because due to charge conservation S^' a =0 
if cr ^ a'. We also notice that, not surprisingly, the first term in the second line can be written 
in a very compact form in terms of the transfer matrix T{q) = str ySipipi, q) ■ ■ ■ Sm,o(pm, q)J , 
a central object in integrable theories. Finally the second term in the last line is irrelevant if 
we only integrate over the S-matrices poles but for relativistic theories we expect the forward 
scattering amplitude to appear when the momenta is integrated over the real axis. 



22 



In Appendix D we consider a direct application of this formula for the SO (oo) sigma 
model. An independent computation of the results therein presented would be useful to put 
the conjecture (I2.38j) in firmer ground. 

In the AdS/CFT setup it would be interesting to consider this expression applied to the 
computation of the exponential corrections to spinning strings [IB! H7] . In the scaling limit the 
transfer matrix becomes the exponential of the algebraic curve quasimomenta [5], El HS? 49, 50J 
and using the techniques in [511 [2H E21 EH [291 G3] one might try to study semi-classical wrapping 
corrections around fairly general classical solutions. 

Finally, we should stress that expression (12.381) is a conjecture for which we have no prove but 
only empirical evidence. It would be interesting to try to directly generalize Luscher arguments 
for many particles states in integrable theories in which factorizability should provide dramatic 
simplifications. 



2.3.3 Perturbative treatment 

In this section we will consider the perturbative small g regime. Since 

e -A» = 2g cos | + 0{g 3 ) 

the leading finite size corrections for generic twists will appear at order g L and at order g 2L for 
twists given by (12. 5p . For example the single magnon results of section [2TT1 are easily seen to 
give 



e(p) - e °o(p) 



— 2JJ ( 2g cos — j cos— (p — O ) , for generic twists 

(2.39) 

—2U \ 2g cos — J , for twists as in (12. 5p 



As explained in section 12.3.11 when we want to consider states with more than one magnon we 
expect two types of contributions. On the one hand, we have a Luscher type contribution of 
the form 

M 

SE Luscher ~ -Uj2 S /3n (2-40) 

n=l 

due to the fact that the energy of the state, as a function of the magnon momenta, changes 
when the state is put in finite volume. On the other hand we obtain extra corrections due to 
the fact that the effective Bethe equation for the magnons are corrected when wrapping effects 
are taken into account, 

se bae ~ 2^ , hn ■ 

n=l dPn 

However, since perturbatively the derivative of the dispersion relation with respect to p carries 
an extra power of g 2 , (I2.40p is sufficient to compute the leading wrapping correction. Notice 
that since the exponential factors of e _/3 °° L start at g L order the prefactors can be computed 
using the g° order positions of the Bethe roots (which we simply denote by u n in what 
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follows) . This is a huge simplification which is probably also present in M = 4 S YM if the most 
fundamental description of the supersymmetric chain bears a resemblance with the Hubbard 
model. 

Generic states can be studied by expanding at small g the expression in Appendix B. For 
concreteness let us focus on 2 magnon states with twists given by (12.51) . Moreover we chose the 
twists as given in [20] which correspond not only to (12.51) but also to e l ^\ L = (— 1) +1 . In this 
case not only the wrapping order is delayed to g 2L but also the BDS effective equations become 
untwisted. We obtain, to leading g 2L order, 



D 2 (AT') + {ui _^ +i +0 ^) (2.4!) 



with 

X'.2 = "J 



(0) _ • L ( 1 \ L U l,2 ~ U 2,l 



u ii2 -i/2J u lj2 -u 2 ,i 



For Konishi like states with opposite momenta p and — p (and thus U\ = —u 2 ) we get the 
remarkably simple expression 

^ = 2 2 ^y^(cos(p) + 3) esc 2 \ cos 2i \ . (2.42) 

In particular, we might consider some " high energy magnons" with p ~ — n in a long spin chain 
to find 

££large momentum ^ + ^2L ^ j 

which is a tiny quantity while for low momentum states with p ~ 

^^ilow momentum ^ _ j 

and wrapping corrections are very large. This is physically intuitive as low momentum states 
probe larger portions of the space cylinder and thus are more sensitive to wrapping interactions. 

It would be very interesting if the computations of [HI [T7j could be generalized to generic 
two magnon states with arbitrary L. If this turns out to be feasible then the study of the 
dependence of the anomalous dimensions on the magnon momenta could be an important 
window into the structure of wrapping effects in M = 4 SYM. 

As mentioned above, to compute ( 12.4111 or (12.421) we only need to know the leading g° 
position of the Bethe roots which are simply given by the Heisenberg chain Bethe equations. 
If «i = —u 2 = u they become simply 

u + i/2\ L u — (—u) + i 



u — i/2j u — (—u) — i 

so that the only effect of the second magnon is to effectively renormalize L to L — 1. Therefore 
the momenta p = 2arctan2ti will be simply given by p = —tc + j^j. E.g. for the Konishi 
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operator p = — ~ and SE = — 2268g 8 which is represented by the mismatch between curves C4 
and C2 in figure [H 

Another example which illustrates the huge variance of the prefactor as function of the 
magnon momenta is the 2 magnon state for some large chain with, e.g., L = 100. For these 
states the wrapping corrections range from the smallest values for the largest momentum states 
with n = 49 for which 8E ~ —7.97 x io~ 300 g 200 to the highest value for the lowest momentum 
states with n — 1 and corresponding 5E = —1.15 x 10 64 g 200 . 



3 Families of long-ranged integrable Hamiltonians 

In the previous sections, we explored wrapping effects in the Hubbard model, as a controlled 
toy model closely related to M = 4, see discussion in the introduction. In that model the long 
range Bethe equations for the magnons are effective equations and the fundamental degrees 
of freedom are electrons whose interactions are ultra local. In this section, we will explore a 
completely different toy model which also has an analytic solution. In this model, the funda- 
mental description is given by a long ranged Hamiltonian where, by construction, wrapping 
interactions are under control. Contrary to the previous model, it does not seem to share many 
features with the known M = 4 spin chain however it is a nice simple model worth looking at. 

The algebraic Bethe ansatz construction is the formalism that diagonalizes transfer matrix 
operators like (11.51) mentioned in the introduction^]. These transfer matrices commute with one 
another for different values of the spectral parameteJ^I. Thus, if we construct a spin chain 
Hamiltonian H from the transfer matrix (usually, by taking derivatives of its logarithm at a 
particular point u*) then, by construction [H, T{u)] = and we immediately obtain a huge 
number of conserved charges and hence quantum integrability. 

For example, let us consider the standard SU(2) spin chain transfer matrix 



f(\) = tr TT — tt^- , (3.1) 




where Voj is the permutation operator between a physical vector space hj and an auxiliary 
space h , both isomorphic to C 2 . The transfer matrix is then an operator acting in the full 
Hilbert space given by the tensor product of L copies h\,...,hi associated with the L spin 
chain sites. The eigenvalues of this SU(2) transfer matrix are given by 

where Uj are denoted by Bethe roots. It is clear from the definition (13. ip that the eigenvalues 
can not have poles as A approaches Uj —i/2. The cancellation of the corresponding residues in 



9 A transfer matrix also appeared in section 12.3.21 although there, the operators being multiplied inside the 
trace were S'-matrices rather than i?-matrices. 

10 The reason for this is that the operators being multiplied in (jl.5|) or in (|3.ip below obey the Yang-Baxter 



triangle relation. 
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(13.21) is indeed guaranteed by the Bethe equations 



Uj + i/2 
i/2 



u ; 



M 



Uj —Uk + l 



Uk 



(3.3) 



which quantize the Bethe roots Uj. Physically, the left hand side of this equation represents 
the free propagation e vpjL of the magnon j as it goes around the spin chain while the r.h.s. 
Ylk^j S(Pj>Pk) describes the scattering of this magnon with all the other magnons. Therefore, 
the Bethe root Uj and the momentum pj of the jth magnon are related by 



u 



1 p 
- cot - 

2 2 



(3.4) 



Having diagonalized T we have automatically diagonalized all Hamiltonians obtained from 
this transfer matrix. For example, the Heisenberg Hamiltonian (I2.6P can be obtained as 



A=0 



so that the spectrum is simply 



E 



1 d 

2i dX 



logT(A) 



A=0 



M 

E 



1/2 



M 



u 2 + 1/4 



5> 

3=1 



sin 



2Pj 



By considering more derivatives of logT at A = we generate other local Hamiltonians with 
longer range. To obtain the spectrum of these Hamiltonians we simple act with more derivatives 
on the logarithm of the eigenvalue (13.21) . In particular, if we take more than L derivatives, 
wrapping interactions, where the range of the Hamiltonian is bigger than the size of the spin 
chain, will appear. Remarkably, for all such models the Bethe equations are just (13. 3p since 
they diagonalize the transfer matrix! 

A particularly interesting hamiltonian is 



H{g) = i log 3^ + h.c. 
Ai s T(0) 



(3.5) 



If we think of g 2 as being an expansion parameter then we have an infinite range Hamiltonian 
where, at each order g 2n in perturbation theory, the interaction range is n. In the notations of 
[56] we have 



4 ^ 



■J +2 J 



+ . . . 



(3.6) 



where T~tj,j+i = 1 — 'Pjj+i- The spectrum of this Hamiltonian is then given by (13.51) where we 



simply replace the operators T(-) by the corresponding eigenvalues T(-) to get 

l m 



M 



E 



Uj-i- g 2 ^ + | 



2 ^ + 47 log 



1 + 



9 



g 2 + i 



n 



U J ~ f - 9 2 Uj-{- g 2 



J=1 U J + 5 



g 2 Uj + | 



+c.c. 
(3.7) 
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The first term comes from the first term in (13.21) . It gives a contribution to the energy of the 
form Yl e (Pj)i that is a sum of the dispersion relations of M individual magnons interacting 
through (13.31) . The dispersion relation, when written in terms of p, yields 



_ 1 / 1 - 2g 2 e^> 



v pI 2 sin § 



2i 



2g2 e +ip/2 g in 



2g sin - + 2g sin p sin - + • 



(2g 2 ) n . np . n p 

H sin — sin - + . . . 

n 2 2 



(3.8) 

The second term in (13. 7p . which comes from the second term in (13. 2p . is identically zero up to 
order g 2L , precisely when wrapping interactions appear! Thus, at order g 2L the energy is given 
by a sum of dispersion relations of the form (13. 8p plus this wrapping term, which entangles 
all M magnons and is not writable as a sum of individual magnon energies. In terms of the 
magnon momentum we have 



M 



M 

L \\ (2e- 2ip i 

3=1 



-i Pj 



) 



c.c. 



+ 0(g 



2L+l\ 



(3-9) 



For example, for 1 and 2 magnons we get respectively 

,2i / / ^ 



E{p) = e{p) + ^ (sm (p+^L\-2 sin ^2p + ^fj ) + 0(g 2L+2 " 



and 



E(p, k) = e(p) + e(k) + g lL s(p) + s(k) 



s(0) 



2s(p + k) +0{g 



J2L+2\ 



(3.10) 



(3.11) 



where s(x) = sin(p + k + Ltc/2 + x). In particular, it is clear that the correction to the energy 
of the two magnon state is not of the form Se(p) + Se(k). This was expected since at order g 2L 
the interaction range covers the entire chain and the notion of asymptotic region where one can 
safely measure the dispersion relation of each magnon is destroyed [91 [55] . 

When we compare the type of corrections (13.101) and (13. lip we get from this long range 
model compared with those we got within the Hubbard model (I2.39P and (12.421) we see that 
while in the latter the coefficient of the order g 2L wrapping correction to the energy exhibited a 
strong L dependence, from factors like cos L (p/2), in the former this coefficient is L independent 
and generically of order 1. If this feature is generic then the weak or strong L dependence of 
the g 2L prefactors to the corrections of 2 magnons states in Af = 4 super Yang-Mills could be 
a sign, respectively, of a fundamental description in terms of a long ranged exact integrable 
Hamiltonian or of a local description in terms of an hidden level of particles in the Beisert- 
Staudacher nested Bethe ansatz equations. In this context, the perturbative computation of (16] 
and [T7], if generalized to a more general setup with 2 magnons and generic L could provide 
very useful hints about the nature of the fundamental integrable structure of M = 4 super 
Yang-Mills. 



27 



3.1 Generalizations and transcendentality 



In the previous section we saw that we could easily generate (long-ranged) integrable Hamilto- 
nians by considering 

1 g 2n d n 

1iH a ^—^ l °sT(\) +h.c. 



H 



n=l 



n\ d\ n 



(3.12) 



A=0 



The spectrum of such Hamiltonians is immediately given by this expression with the operator 
T(A) replaced by the corresponding eigenvalue ( 13.2H . At order g 2n these Hamiltonians are 
local with interactions of range n. If we truncate the expansion at a given order m by setting 
«n>m = 0, then for chains of length larger than m we have no wrapping interactions and the 
energy is simply given by a sum of dispersion relations ^ . e{pj), wit 

, . l^a n g 2 U 1 1 



e it 



2% ^ 

n=l 



n 



(u-i/2) n (u + i/2) n 



(3.13) 



If, on the other hand, we consider an infinie sum or, alternatively, if we truncate the expansion 
in g 2 at an order m > L, the spectrum (starting at wrapping order g 2L ) will no longer be a 
sum of individual dispersion relations. In particular, precisely at order g 2L we obtain 



M 



E 



2L 



M 



n 



Uj — 3i/2 
i/2 



c.c. 



+ 0(g 2L+2 ) 



(3.14) 



where the second term takes into account the wrapping interactions. It is interesting to notice 
that for all these families of long ranged Hamiltonians the expression for the wrapping inter- 
actions is quite simple and absolutely universal. The example f !3.5j) we considered corresponds 
to a n — 1 but as we see any choice of a n will lead to a solvable problem. A particularly funny 
example would be 



H 



1 oo 

-Yc n 

2i^ 



9 



2n 



d 



■2 n 



n=l 



(2n- l)\dX 



2n 



logT(A) 



+ h.c. 



A=0 



with C n the Catalan numbers, for which we find the curious expression 



eiu 



X+(u) X-{u) 



u = 
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for the dispersion relation. As a function of the Bethe roots this is precisely the dispersion 
relation appearing in the BDS equations jl] and even in the full AdS/CFT Bethe equations [8]. 
However, unfortunately, as a function of the magnon momenta this is not the same as (12. lOj) 
because the relation between u and p in our toy models is always of the form (13 .4[) instead 
of (12.121) . In other words, although the desired dispersion relation can easily be obtained, the 
Bethe roots are always quantized via the usual Heisenberg spin chain Bethe equations. 



11 If the a n are not real - which from the definition 13. 121 is a perfectly reasonable possibility - we should take 
the real part of this expression. 
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We can generalize this construction of long ranged Hamiltonians for other symmetry groups 
as well. For example, for SU (N) the transfer matrix in the fundamental representation takes 
exactly the same form (13.11) as for the SU(2) chains except that the operators being multiplied 
now live in hj ® h where both h Q and hj are copies of C . Thus, we can still consider 
Hamiltonians of the form (13. 5p with a g 2 perturbative expansion (13.61) . As before, they will be 
long ranged Hamiltonians where the Hamiltonian range at order g 2n is n. The only thing that 
changes is the expression for the eigenvalue (I3.2p . For an SU(N) spin chain we have N — 1 
types of roots and K n Bethe roots of the type n = 1, . . . , N — 1. The transfer matrix eigenvalue 
reads^ 



Ki \ (1) _ • /9 / , \ L N—l K n \ M n +2- K n +i x _ (n+l) n _l . 



/ 



We see again that when we compute the first few local charges expanding logT(u) around u = 
only the first term gives a non-vanishing contribution. As before we can consider Hamiltonians 
of the form (13.121) . The dispersion relation for the SU(N) magnons is exactly the same as for 
SU{2) while the leading wrapping correction to the spectrum can be computed as before to 
yield the generalization of (13.141) . 



Ki 2L 
3=1 



N-l K n (n) _ n +2 ■ K n +i (n+l) _ n-1 • 



n=l 3=1 "j 2 3=1 j 2 



where again P = 4 Y] . log is the state total momentum. 

& i & uf>-i/2 

All these considerations can be trivially generalized both to non-compact spin chains and 
to supersymmetric ones. For example, for the SL(2) spin chain we have [53, EH] 

r -«M=ni:£+$ + |j(sTx) J 8(A-«. J + ( ^o(A , -i+a?ii) (3 - 15) 

with the Bethe equations, following from canceling the poles in this expression, reading 

u j + i/2 \ L ™ u 3 -u k -i 
u 3~ % l 2 ) t£ u i~ Uk + l 

which differ from ( 13.31) by a simple sign in the r.h.s. Again, by expanding the log of the transfer 
matrix around A = we see that only the first term contributes until the L'th derivative is 
taken. Thus if we consider an Hamiltonian of the form (13. 12[) we will have, up to order g 2L , 
the energy as a sum of the same dispersion relations (13.131) . In particular if the constants a n 
are algebraic numbers then so will be . e(uj) when truncated to order g 2L because clearly the 
solutions to (I3.16P are also algebraic (complex) numbers. 
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The SU(N) Bethe equations can be immediatly obtained by canceling the apparent A poles in this expres- 
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However, precisely at order g 2L the second term in (I3.15P starts contributing and we find 



M 




;2L 



oo M 



E = Y J <u ] ) + 



9 



(3.17) 





— c.c. 



n=l j=l 



This new wrapping term differs from the one computed for the compact groups SU (N) by the 
fact that it is given by an infinite sum of terms. Thus even if Uj and a n are perfectly algebraic 
numbers the energy of this state will only be algebraic up to order g 2L , when this infinite sum 
will give a transcendental contribution! 

Let us consider a few examples. We chose ai = 1 for simplicity. For L = 4 the one magnon 
state with momentum 27r/4 will be corrected to 



In (11.81) in the introduction we listed a couple of additional examples. 

Although, this model is not immediately related to the (non-compact sector of) AdS / CFT 
Bethe equations which are much more complicated than (I3.16p . it is still interesting to see 
that transcendentality very naturally appears due to the non-compact nature of the transfer 
matrix. In particular, if an extra level of hidden degrees of freedom is to be discovered then 
the transcendental numbers present in the dressing factor could be an important hint. A more 
fundamental PSU(2, 2|4) symmetric transfer matrix in the field strength representation would 
also be given by an infinite sum of terms since the representation is infinite dimensional. It 
would be spectacular if a relatively simple extended transfer matrix with some extra degrees of 
freedom included and only simple algebraic expressions could lead to the intricate structure of 
the full dressing factor where transcendental numbers abound. Probably the correct place to 
try to reverse engineer and find this extra level of hidden particles is the transfer matrix rather 
than the Bethe equations. 
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13 For two magnon with opposite momenta the sl(2) equations are trivially solved exactly as explained in the 
end of section l2.3.3l for the su(2) chain. For the non-compact chain the effect of the second magnon is simply to 
renormalize L — > L + 1 instead of L — > L — 1 as we had for the SU(2) chain. Thus we obtain pi = —p2 = j^j 
for the two magnon state with opposite symmetric momenta. 




2C(3) + 2C(5)) + 0G ? 12 ) 




we 
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Appendix A Wave function of the Hubbard magnon 

In this section we shall study states made out of two fundamental particles o and |, 

|*) = 5>(n,n') (< : 4, f0 ) | T • - T> ■ 

n,n' 

Acting on this state with the Hamiltonian ( 12. ip we can find the form of the wave function 
if)(n,n') so that the state is an eigenstate, 

H\V) = E\V) . 

This gives the equation 

Eip(n, ri) = -t [e^ipin + l,ri) + e'^{n - 1, ri) + e^°ip(n, ri + 1) + e^ipin, ri - 1)] 
(E + U) ip(n, ri) = -t [e^ipin + l,n) + e -i ^(n - 1, ri) + j*°ij)(n, n + 1) + e~ i<t>0 ip(n, n - 1)] 

with ri ^ n. A plane wave superposition 

^(n,ri) = e -ih n -^° n ' (Ae iqn+iq ' n ' + Be iq ' n+iqn '^ , n < ri (A.l) 
i/;(n,ri) = e -%»-^«' (Ce****"' + De iq ' n+iqn '^ , n > ri (A.2) 
solves the first equation and yields 

E = — 2t cos(g) — 2t cos(g') . 
Continuity of the wave function at n = ri gives 

A+B = C+D 

and the second equation reduces to 

A (2g sin(g) - 2g sin(g') - i) = C (2g sin(g) - 2# sin(g')) + iB . (A.3) 

In an infinitely large volume we might look for bound states with exponentially decaying wave 
functions. These are only possible if two of the exponentials in flA.lf) and (1A.2|) disappear. 
Parametrizing the momenta as q = p/2 — i/3 and q' = p/2 + ij3 we see that we need B = C = 
and A = D ^ 0. Then equation (1A.3[) fixes 

sinh /3 = sinh (3^ = (A. 4) 

Ag cos | 

in which case 

tfj( n ,ri) = e-^-^n' e ^(n+n') e -/3\n-n'\ 



and the energy (1A.3I) reads 



E = £oo (p) = —At cos - cosh /?oo = — Uy 1 + 16c/ 2 cos 2 ^ . 



31 



Comparing with the energy (1A.3j) of two particles with momentum p/2 we conclude that the 



magnon has a relative binding energy 

cosh (3 — 1 . 

On the other hand, at finite volume we impose periodicity 

i[)(n + L, n) = ip(n, n) , n < n < n + L 

and 

ip(n, n + L) = ip{n, n) , n < n < n + L 
of the wave function to obtain 

Ce i{q-<f>l)L = A ? De iW-^)L = B 

and 

Ae i{ q '-4>o)L = c ^ Be i(g-ML = D 
This immediately gives total momentum quantization 

e i(q+q' ' -4>i-<t>o)L _ ^ 

and the relation 

C ~ 1 - e i(q~<t>o)L 

From ( 1A.3I) we obtain 

l[q - h)L _ 2^sin(g)-2^sin(gO + z ;: e J ( ( ;.^ 



2(yfsin(g) — 2gsm.(q') — i 

which can be rewritten as 



^ (sin(g) _ sin(g')) (l + e^i'^) + yj g 2 (sin(g) - sin(g')) 2 (l - e i{ ^-*° )L ) 2 - e^x-^ 

2g sin(g) — 2g sin(g') — % 

This is precisely what one obtains from the Lieb-Wu equations ( 12. 4ft for the Bethe roots q,q',u 
eliminating the auxiliar variable u using its Bethe equation. Moreover, in the case of identical 
twists this becomes the simple S-matrix 



e i(q-<f>o)L 



2g sin(g) — 2g sin(g') + i 



2gsin(g) — 2gsin(q') — i 

If again we look for solutions in the form q = p/2 — i(3 and q' = p/2 + i/3 so that 

p 

E = e(p) = —At cos - cosh (5 , 

we obtain 

jto-io-'H) 1 ' = 1 (A.5) 
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and 

a P • i o sinh(/3L) 
Ag cos - sinn p — 



2 M cosh(pX) -cosL (§ -</>„) ' 
Notice that the result is symmetric under the exchange of twists since (1A.5I) implies 

cosL (| - = cosL (~ - 0|) = cos ^ (0 : - (f> ) . 

In the case of equal twists the expression reduces to 

p . sinh(/3L) 
Ag cos - sinn p - 



2 ^ cosh(/3L) ± 1 

and for BDS twists obeying e tL ^°~^ = —1 gives 

p 

4g cos — sinh (3 = — tanh (pX) . 

In any case, for large L the r.h.s. can be replaced by —1 and we recover the infinite volume 
result (IA.4I) . 



Appendix B Dressed BDS BAE with twists 

This appendix complements the computations and results of sections (12.3.11) and (12.3.31) . As we 
see from (12.341) . (12.361) and (I2.32f) a crucial quantity we need to compute the leading wrapping 
corrections is A n defined in (12.311) . To compute A n we focus on the equation for a single 
constituent of the Takahashi triplet, say on equation (12.231) for q n . Expanding this equation 
using (I2T311 and fl2T33|) yields 

A n ~ M°) = le i E f+^~'t>o)L TT U m-U m +l = r-^ L) (R1) 

™ J. J. 7 ,00 ? ,00 \ J 1 \ J 

m+n m m 

to leading order. However, as seen from the above mentioned equations (12.341) . (12. 36ft and (I2.32H 
what we really need is the imaginary part of A n . Due to the corrected BDS equations (I2.34[) 
the imaginary part of A„ is of order e~@ nL for generic twists and of order e~ 3/3ni when (12.51) . 
Indeed , the imaginary part of Ai ' can be simplified using the BDS corrected equations (I2.34[) 
to give 



3 



(Af) ~ -iAf e «(*-*H cos^(0 o - : ) (B.2) 



which indeed vanishes when (12. 5p . Thus, in this case we need to expand the equations for q n 
further, 

/ ■> A (0) A (0) \ 
A n ~ AW - i (AW) 2 + V %^ . (B.3) 
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To find the corrected BDS equations we multiply the equations for the Takahashi triplet to get 



,00 „,oo • „,oo „,oo 



where 



ji(pn-<j>i-<t>o)L U n U m \ _ U n U m \ U n u m 1 r m,n r m,n+M ^ 

m ^ n U n ~ U m ' % m=1 U n ~ U m ' 1 U n ~ u m + * r n,m r n,m+M 



^ u n~ 2gsmq m - i/2 
n ' m u n -2gsmq m + i/2' " ' 



Notice that so far no approximation whatsoever was done and (1B.4|) is an exact relation. Ex- 



panding now the r.h.s. of (1B.4|) using (I2.3ip and (I2.33P we find the corrected (twisted) BDS 
equations (I2.34p 

Again, when the twists are chosen as in (12. 5p we need to be more careful and expand our 
expressions further. By expanding the product of the equations for the Takahashi triplet to 
order e~ 2/3L we get the same (I2.34p with 

^(^n) ( 1 , 9 Vn . 2 \ 2[U n W m )A n A m . . . 

tan 1 + — — — — (n m) (B.6) 



(m„ - M m ) 2 + 1 V«„ 2 u n -u m J ((u n - u m ) 2 + 1) 

instead of (12351) . 

Finally, in (12.361) the energy is corrected because the momenta change (the second term) and 
because the functional dependence on the momenta changes (the first term). For completeness 
let us mention why the correction due to the change in functional dependence is so simple. This 
is simply because 

- to cos ^ cosh (A^) +5p n ) = -U ^-40cos^sinh/3 M (p~)J 5(3 n = -U5f3 n (B.7) 
since — Ag cos -g- sinh ^(p^) = 1. 



Appendix C Luscher with magnon-magnon S-matrix 

at strong coupling 

In this appendix we expand the results of |27j^l concerning the BDS Luscher term to the next 
order in 1/g. The computation in question is the /i-term 

= «** I'M -'<**)] z^y + - 

14 The notation in this paper differs from ours by pthere = Phere + ir, gthere — V^ghere- Moreover, we have 
an extra factor of U multiplying the magnon dispersion relation. This shift in p introduces an annoying (— 1) L 
factor. In order to avoid carrying it around we consider L even in this appendix. 
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where 



u (p) = 2 tan § J 1 + 16 9 2 cos2 | 



and k* is defined by the pole condition 

u(p) - u(k*) + 2 = 0. 

For large g, this gives 

k = 7r H sec — tan -sec — (23 cos p — 49) sec — h... 

2g 2 4g 2 2 2 384g 3 v F ' 2 

which yields 

xf\ sof U P, iU , P 2P, \ L P(, i P P, 

oelp) = yi sec — — - tan - sec — h • • • exp sec - 1 tan - sec — h • • • 

K } V 9 2 g* 2 2 J P 2g 2 \ 2g 2 2 

Appendix D Many particle Luscher formula and 

SO (oo) sigma model 

In this appendix we apply our conjectured expression (12.381) to the study of finite size corrections 
to the energy of many particle states in the SO (oo) sigma model. The motivation for such 
computation is twofold. On the one hand, it is a very interesting model where the single 
particle corrections are very well understood [59j [311 [32]. On the other hand, large N sigma 
models often allow for the computations to be done in a vast panoply of ways and in particular 
it might be simple to cross check our conjecture through some saddle point computations. 

To apply Q2.38P we need to compute £%l'£(pi, q) S^fcfai, q) ■ ■ ■ s IZ'Im^ l) which is P re " 
cisely the transfer matrix T(9; {0j}) = tr ^Si,o(0i — 9) . . . Sm,o(Qm — Q)\ if we parametrize 
Pj = msmh9j and q = msinh^. In the large N limit the S'-matrix tends to unity but since 
we will sum over a large number N of such terms, the sum we need to compute converges to a 
nice quantity, namely, using the normalizations of [32], 

M „ . 

S a SM, q)SZ%(pz, q)... S£,t>, q) - # . . • - £ "g, 

j=l 3) 

where we sum over the index of the virtual fundamental particle bj. There are no bound 
states in this theory and therefore the /z-term is zero. Thus, we will consider only F-term 
like contributions. Finally we have dq(e'(q) — e'(pj)) = md6 sinh(# — 9j)/ sinh^ and e tqL = 
e imLsmhe . Will all ingredients at hand it is then straightforward to apply f!2.38j) to find 



for a single particle state and 



M r M 

6E(9 1 , 9 M ) = + E 5e ^t(0i, 02) (D.2) 

j=l ■? i<j 
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for a many particle state. Here 

S J1 = f d6e- mLcoshe = 2K (mL), 
m J 

in agreement with [59], and 

S Jj =m f dQe-mLco^e ( cosh(6 - 6j) cosh(fl - OJ) \ 

int J Vcosh(^) cosh(fl - 6j) cosh(^) cosh(# - Oi)) 1 ' 



These expressions have a clear physical meaning. As for (ID.lj) we get the expected dependence 



on the momenta since we know that for a single excitation the form of the dispersion relation 
will change solely due to the change in the particle mass when the system is put in finite 
volume and thus we expect 8e = 5\J m 2 + p 2 = mdm / a/ m 2 + p 2 = 5m/ cosh(6 l ) precisely as 
found. For the many particle particle state we get expression flD.2j) where we see that the energy 
of the state is corrected due to the mass shift of each particle and by an extra interaction term 
flD.3j) . Probably this interaction correction to the energy of the many particle state could 



be independently determined through some sort of saddle point computation. This would be 
extremely interesting as it could help to partially confirm (or correct) (12.381) . 
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